Parametric instabilities in magnetized multicomponent plasmas 
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This paper investigates the excitation of various natural modes in a magnetized bi-ion or dusty 
plasma. The excitation is provided by parametrically pumping the magnetic field. Here two ion-like 
species are allowed to be fully mobile. This generalizes our previous work where the second heavy 
species was taken to be stationary. Their collection of charge from the background neutral plasma 
modifies the dispersion properties of the pump and excited waves. The introduction of an extra 
mobile species adds extra modes to both these types of waves. We firstly investigate the pump 
wave in detail, in the case where the background magnetic field is perpendicular to the direction 
of propagation of the pump wave. Then we derive the dispersion equation relating the pump to 
the excited wave for modes propagating parallel to the background magnetic field. It is found that 
there are a total of twelve resonant interactions allowed, whose various growth rates are calculated 
and discussed. 

PACS numbers: 52.35.Bj, 52.35.Mw, 52.30.Ex, 52.25.Vy 



I. INTRODUCTION 

The basic natural modes of magnetized plasmas such 
as those that occur in molecular clouds, cometary plas- 
mas and stellar atmospheres are of great interest. When 
the frequencies are low then a class of linear waves, re- 
ferred to collectively as Alfven waves, are known to exist 
and are of importance to the understanding of many ba- 
sic plasma phenomena P, Q ■ The linear approximation 
to these waves breaks down at large amplitudes where 
nonlinear effects become important in their propagation. 
One such large amplitude wave is a magnetoacoustic 
wave, which modifies the background magnetic field in an 
oscillatory fashion, and so can be considered as a pump 
wave that drives other waves nonlinearly. Such large- 
amplitude pump waves may occur in conditions such as 
seen in solar and space plasmas. For example, solar shock 
waves can set up large amplitude standing magnetoacous- 
tic waves in coronal loops or magnetic flux tubes 0,01 HI- 

By considering perturbations of this large amplitude 
pump, we can investigate the possibility of exciting nat- 
ural modes of the system, such as Alfven waves, due to 
a resonant interaction between the pump wave and the 
natural waves. In the single ion species case this basic 
phenomenon was predicted by pi H , and subsequently 
pursued by several authors |8|, |3, [l0| . For a pump mag- 
netic field parallel to the background magnetic field in a 
single ion species plasma, the excited waves are Alfven 
waves travelling in opposite directions along the magnetic 
field 0,0- If ion-cyclotron effects are included, the ex- 
cited waves are the fast and slow (ion-cyclotron) Alfven 
waves travelling in opposite directions 8]. More general 
behavior is allowed in the case where the excited waves 
are permitted to travel obliquely to the magnetic field 

m 
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The parametric excitation of waves in a dusty magne- 
tized plasma has been investigated in Ref . [l2j , but in the 
approximation where the dust is taken to be immobile. 
In that case the oppositely travelling pair of waves are 
modified, with the presence of a cutoff frequency in the 
fast Alfven wave. Furthermore, an interaction between 
a pair of fast waves or a pair of slow waves was found 
to exist. The corresponding growth rates of the slow- fast 
and fast-fast pairs were maximized as a function of the 
dust concentration. 

Many space and laboratory plasmas are multicompo- 
nent, i.e., contain multiple ion species. It is therefore 
of great interest to investigate the effects of the addi- 
tional species on the linear and nonlinear properties of 
the waves in the plasma. For instance, it is often the 
case that in the presence of an additional ion compo- 
nent an extra mode is excited, or forbidden regions of 
frequency are introduced. One area of study of multi- 
ion plasmas is the bi-ion plasma. The inclusion of one 
extra ion species often captures the basic information of 
several extra species. The bi-ion plasma has particular 
importance in plasma fusion, laboratory plasmas and in 
astrophysical environments, with the secondary ion usu- 
ally positively charged. 

A dusty plasma adds another level of interest to the 
topic of multicomponent plasmas, due to the dust prop- 
erties. Dust is an additional impurity of large mass and 
often of negative charge. Under the simplest approxima- 
tion, all the dust grains may be considered to have the 
same mass and equilibrium charge, and so are equivalent 
to a second ion species. This is the case considered in this 
paper. It is known that the inclusion of dust in a plasma 
may introduce cutoff frequencies into the basic Alfven 
waves, and introduce a low frequency mode, whose na- 
ture is different from that in a bi-i on p lasma due to the 
dust grain's extremely high mass [T^, [Tj|- Qualitative 
differences to the bi-ion case arise when charge perturba- 
tions of the grains are included, leading to an additional 
damping mechanism or when a spectrum of dust 
grain sizes and charges is allowed for [l(j- Dust is found, 
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in varying amounts, in many astrophysical and space en- 
vironments such as molecular clouds and the rings of Sat- 
urn. Dusty plasmas have also been studied closely ex- 
perimentally, since the heavy grain mass introduces low 
frequency effects that may be studied in real time. 

An immediate consequence of the presence of an ad- 
ditional ion species is its modification of the background 
free-electron number density. This influences the propa- 
gation of plasma and electromagnetic waves [lj, Ull 03 • 
If in the case of a dusty plasma we fix the grain charge, 
then the grain becomes entirely defined in terms of its 
mass and charge, and as such acts just like another ion 
in the plasma (albeit of negative charge). The state of 
charge neutrality may be written as: 



eZ\n\ + eZ 2 n 2 = 0. 



(1) 



Here 7i e ,i,2 is the number density of plasma electrons 
(with the charge — e), and the two ion species (of signed 
charge Z\ and Z 2 ), respectively. For laboratory dusty 
plasmas, the grain charge is negative (i.e., Z 2 < 0) and 
large (\Z 2 \ ~ 10 2 — 10 3 ), so that an appreciable propor- 
tion of the negative charge in the plasma may reside on 
the dust particles. For astrophysical dusty plasmas, \Z 2 \ 
may only be of the order of unity, which is often the case 
for a canonical bi-ion plasma. In environments such as 
the interstellar medium, where the dust grain is in an 
electron-proton plasma with little ultra-violet radiation 
present, the dust grains acquire a negative charge. On 
the other hand, exposure to ultraviolet light from nearby 
stars can cause ionization of the grains, leaving a residual 
positive charge. 

In this paper, we investigate the propagation of plane 
hydromagnetic waves (Alfven waves and magnetoacous- 
tic waves modified by the presence of a secondary ion 
or dust species in a bi-ion plasma), propagating paral- 
lel to the pumped magnetic field of a large amplitude 
magnetoacoustic wave. We generalize Ref. [lj] where the 
second heavy space was assumed immobile, to the case 
where both ions are fully mobile. A further generaliza- 
tion is the inclusion of pressure. The background mag- 
netic field and plasma density are taken to be uniform, at 
frequencies well below the electron plasma and cyclotron 
frequencies. First, we find the dispersion equation of the 
pump wave, and then concentrate on pump waves of large 
wavelength. We then obtain a coupled pair of equations 
of motion governing the perturbed plasma. Resonant 
interactions are sought, and growth rates of the waves 
parametrically excited by the pump are calculated and 
discussed. 



II. MULTIFLUID MODEL 

The most general set of equations used to describe 
two mobile ion-like species, plus electrons, is a three-fluid 
model. In this picture we employ three momentum equa- 
tions for the electrons and the two ion species, where we 
include both the ions' inertia terms, while ignoring the 



electron inertia. This is valid if the frequencies of interest 
are well below the electron cyclotron frequency. In addi- 
tion we use Maxwell's equations and two mass continuity 
equations for each of the ion species, ignoring Maxwell's 
displacement current. 

The primary species (assumed positively charged) shall 
be denoted with a "1" subscript, and the secondary 
species (either positively or negatively charged) shall be 
denoted with a "2" subscript. In this notation we may 
write the current density as 



J = Z\en\V\ + Z267J2V2 — en e v e 



(2) 



where the two ion species and electrons have velocities 
denoted by Vi, V2, and v e , respectively. We employ the 
parameters <5i ~ n e jZ\n\ and S 2 = n e /Z 2 n 2l which mea- 
sure the distribution of charge in the plasma amongst 
the ions. Employing the total charge neutrality condi- 
tion, given in Eq. Q , we may write this as 



1 1 

T + T = 1 ' 

Ol o 2 



(3) 



in terms of these parameters. Note that in the limit of a 
single (primary) ion species, we have 5i — > 1 and |5 2 | — > 
oo. 

Ignoring collisions, but including the effects of pres- 
sure, the starting equations for the velocities, electric 
and magnetic fields and each number density are the mo- 
mentum equations, the two ion continuity equations, and 
Ampere's law neglecting the displacement current: 



mini 



m 2 n 2 



dvi 

~dT 
dv 2 

~dT 


dni 

~dT 

dn 2 
~dt 



Here, Pi t 2,e are the two heavy species and electron pres- 
sures, and mi^ are the heavy masses. The magnetic field 
B includes the background magnetic field Bo- 

Finally, by assuming either an isothermal or an adia- 
batic equation of state, we have 



-Vpi + Zieni (E + vi 


xB), 


(4) 


-Vp 2 + Z 2 en 2 (E + v 2 


xB), 


(5) 


—Vp e — en e (E + v e x 


B), 


(6) 


-V-(mvi), 




(7) 


-V • (n 2 v 2 ), 




(8) 


li e(Zimvi + Z 2 n 2 v 2 


- n e v e ). 


(9) 



Vp Q = U a m a Vn a 



(10) 



for each species, where U a are the individual sound 
speeds. Though it is not imperative, we shall assume 
isothermal changes, which permits us to write down U a in 
terms of the plasma temperatures, i.e., U\ — kBT a /m a 
is the square of each thermal speed, where T a is each 
temperature and fee is Boltzmann's constant. 

At this point we choose to eliminate the electron vari- 
ables from Eqs. 10}-©, and use Faraday's law to elimi- 
nate the electric field E. There are then two choices for 
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the way to proceed. We can add and subtract the mo- 
mentum equations and deal with a total fluid velocity 
and the current density J. Instead we shall employ the 
most direct method, which is to deal with the ion veloci- 
ties Vi and V2 separately, since our equations will exhibit 
the most symmetry this way. By employing the charge 
neutrality condition to all orders, the following system of 
nonlinear partial differential equations is found: 



d P 2 V7 I \ 

^ = -V-(p 2 v 2 ), 

OB B dvi 

— = Vx(v 1 xB)--Vx — 

= Vx(v 2 xB)-^Vx^, 
1I2 at 



(11) 
(12) 

(13) 



-(*iE7? + a?)Vpi-alV/>a 

^(vi - v 2 ) x B + — (V x B) x B, (14) 

i>0 MO 



S 2P 2^ = -(S2U! 



al)Vp 2 - afVpx 



+ ^(v 2 - vi) x B + — (V x B) x B. (15) 

In Eqs. l(TT |l -H15 |) pi ;2 = "ii, 2 ni, 2 are the densities of 
each massive component of the plasma and f2i j2 = 
Zi t 2eBo/mi t 2 are the corresponding (signed) cyclotron 
frequencies, with Bo = |Bo|. Also a\ — Z\U%m e /mi and 
oi\ = Z2U^m e /m2 are pseudo-squared thermal speeds as- 
sociated with the electron pressure (if Z 2 < we allow a 2 
to be imaginary) . The presence of both Vpi and Vp 2 in 
both equations of motion, is a consequence of the elimi- 
nation of the electron variables from the equations. To be 
more explicit, the electron density fluctuations generate 
fluctuations in both the ion densities through the charge 
neutrality condition, which produces this coupling in the 
momentum equations. The second term on the right 
hand side of Eq. (|13JI is the Hall term; it is important 
when the wave frequency is comparable to either of the 
cyclotron frequencies. Also, we see that the two species 
are strongly coupled through the momentum equations, 
via the third term on the right hand side in each. This is 
an additional Hall-type term associated with the relative 
motion of the two species. Note Eqs. <|11[) - (|14[1 reduce, in 
the case of a single species plasma (Si — > 1, \S 2 \ — > oo, 
assuming local charge neutralit y is maintained), to the 
equations used in Ref. [H HI lU where non- 

linear Alfven waves were investigated. They are then 
known as the (collisionless) Hall-MHD (magnetohydro- 
dynamics) equations. 

The neglect of the displacement current in Maxwell's 
equation is justified when the electron current is much 
greater than the displacement current, which leads to 
the conditions 



J A2 



in 2 



(16) 



where u> is the wave frequency and vai , vai are the Alfven 
speeds associated with each ion. The charge neutrality 
condition implies that the wave frequencies are restricted 
to the regime u) <C uj pe , which may be restated as 



1*1 



» 



jOi, 



J A1,2 



(17) 



There exists a wide range of physical environments for 
which all of the above mentioned conditions are met. 



III. MAGNETOACOUSTIC PUMP WAVE 

Suppose there is a constant background magnetic field 
in the z-direction, given by Bqz. We now periodically 
pump the field, with a periodic modulation: 



Bq[1 + £&cos(w i)]z, 



(18) 



where ujq is the pump frequency. Here e is a constant 
dimensionless quantity which determines the amplitude 
of the pump. More precisely, e is an expansion param- 
eter, which permits us to keep track of our terms, by 
matching powers of e. The parameter b on the other 
hand is a dimensionless quantity that wc include to cap- 
ture any necessary frequency information, i.e., b = fe(^o), 
such that the average of b over ujq is 0(1). By speci- 
fying a particular choice of normalization condition, we 
may solve for b. We later in fact impose the condition 
that the energy density in the pump system is a constant 
over ljo, and calculate the resulting b. Since wc are in- 
terested in a large amplitude pump wave, s will typically 
be O(10- 2 - 10- 1 ). 

However, since this field has no spatial dependence, it 
does not satisfy our wave equations. Hence this is an 
approximation to a wave with a large wavelength. Wc 
must therefore modify the pump with an envelope of 
some wavelength. Since the pump wave magnetic field 
points in the z-direction, the wavevector should point 
perpendicular to this to satisfy V • B = 0. For a planar 
geometry plasma, we choose the axes such that the wave 
varies in the x-direction, with wavenumber ko- By denot- 
ing the pump magnetic field (and all subsequent pump 
fields) with a "0" superscript, we have 



B (0) =B [1 + ebcos(k x) cos(u i)]z- 



(19) 



In a cylindrical plasma, with r denoting the radial dis- 
tance from the cylindrical axis, we have 



B (0) = B [l + ebJ {k r) cos(w t)]z, 



(20) 



where Jo is the Oth order Bessel function of the first kind. 
Note that b = b(k , w ) in Eqs. (O and (|2U|) . 

In the absence of a wave, we shall suppose that the 
plasma is stationary, with ion densities poi,2 and charge 
ratios of <5qi j2 . The effect of pumping will be to modify 
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the velocity, density and charge imbalance to order s, so 
that 



£Vi, 



ev 2 , 



(o) 

Pi = Poi 



epi, 

#i = S 01 + eSi, 



(o) 

P2 = P02 
X<f>) - . 



" E P2i 

eS 2 . 



(21) 
(22) 
(23) 



These quantities define the pump wave, and we now pro- 
ceed to solve the equations (|11[1 - H15|1 to order e (i.e., the 
linear solution). 

By inspecting the form of the magnetic field in Eq. I|19|) 
[or l(2*D)l]. and Eqs. (T3J| and JT5J, we note that the velocity 
cannot have a z-component. Hence, the velocities of the 
ion species must be of the form 



Vi = (v x i,v yl ,0) , v 2 = (v x2 ,v V 2,0) 



(24) 



for planar waves. From the resulting pair of equations 
we ascertain the form of each velocity component and 
density: 



v x i,2 = A xl>2 sm(k Q x) sm(ui t), 
Vyi.2 = sin(fc a;) cos(w £); 

Pi, 2 = Ri,2 cos{k x) cos(w i), 



(25) 
(26) 
(27) 



where Aci,2j an d -Ri,2 are amplitudes independent 

of space and time, to be determined. Then, from the 
continuity equations, we have 



-f»i — An — A x \, 

UJ 



R.2 = P02—A x2 . (28) 

UJQ 



For the cylindrical plasma case, we make the substi- 
tutions: x — ► r, ?/ — > 0, cos — * Jo and sin — > J\, where 
4> is the azimuthal angle and J\ is the 1st order Bessel 
function of the first kind. Then it is found that all the re- 
lationships pi )l -P% jl still hold. It follows that the modes 
in the planar and cylindrical plasmas satisfy the same 
dispersion equation. 

Upon substitution of Eqs. J23)-(|2EJ into Eqs. fTTJl - 
(|15f) . written to first order in e, the following dispersion 
equation of the pump wave is obtained: 



2/2 

uj {uj - 
Here, we have defined 



= Wk 2 (uj 2 -X-Yk 2 ) 



UJ C = Qi/Sqs + Q2/S01, 
W = v\ 1 {l + 5l l l3 l + § Ql B l )/8l 1 
+ v 2 A2 {l + 5 2 2 p 2 + 8 m B 2 )/5 2 2l 
X = u c {v 2 A1 VL 2 (l +0i + SoiB^/Sq! 

Y = v 2 A1 v 2 A2 (p 2 (l + 5oiB 1 )/6 2 01 
+ /3i(l + 5 02 B 2 )/S 2 2 + fcfa)/W, 



(29) 

(30) 
(31) 
(32) 
(33) 



where \lj c \ is a hybrid cutoff frequency as fc ~~ * 0. 
Also we have introduced j3\ := Uf/v A1 , (3 2 := U 2 /v A2 , 
B\ := a\/v 2 Al and B 2 := a 2 /v A2 which allow us to 
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FIG. 1: The normalized pump frequency wo/fii versus nor- 
malized pump wavenumber v A iko/£li with S01 = 1.1 and 
/ f^i = 0.1. (a) The warm plasma, with f3i = 0.8, Bi = 1.5 
and Ua/Ui = 1/2. (b) The cold plasma, with f3i = Bi = 
U 2 = 0. 



specify the cold/warm plasma regimes, with B\ and 
B 2 related through: &q\B\ — Sq 2 B 2 . Note that as 
5qi — * 1 , I <5o2 1 — > 00 (i.e., a single species) we are left 
with the nondispersive relation 



CJ 



+ v 2 A1 k 



(34) 



where c 2 — U 2 + a 2 is the combined sound speed in 
this limit. This is the familiar fast magnetoacoustic wave 
characteristic for k _L Bo- However, in the presence of 
a secondary ion species the fast magnetoacoustic wave 
gains an additional mode and the relationship is disper- 
sive, see Fig.^ 

It is easy to show that the parameters X and W are 
always positive, while Y = for the cold plasma and is 
positive if any of fi\ , f3 2 , B\ , B 2 are nonzero. Although 
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Fig. ^ is for a plasma in which the secondary ion species 
is positive, the same basic qualitative features are present 
when the second ion species is negative. For the warm 
plasma the effect of the Yk^ term in the dispersion re- 
lation is to cause the lower branch to increase without 
bound as fco — » oo, as indicated in Fig.^a). Physically, 
it can be thought of as the fast magnetoacoustic mode 
being converted into an acoustic mode for large fco due to 
the inclusion of pressure. However, if the plasma is cold 
(i.e., Y = 0) then the lower branch experiences a reso- 
nance as fco — > 00 j as indicated in Fig.^b). In this limit 
X becomes the square of a resonance frequency, given by 



X 



2 



(35) 



which gives rise to a forbidden frequency region between 
|uv| and |w c |, since \u> r \ < \u> c \ whenever Sqi 7^ 1. 

The ratio of the cc-velocity amplitudes A x \^ to the 
magnetic field parameter b are of particular importance 
to our later discussion, so we briefly discuss them here. 
In terms of the frequency and wavenumber, we find, 



(36) 



where 



A 2 



A x i _ 


u> 


A x2 


u> 




b 




b 


k A 2 ' 




1 


1 n 2 


U 2 k 2 


+ Q,lLU c — 


UJ 2 


Sot 


$02 fil 


n 2 k 2 

u 2 K 


+ ri2UJc — 


u 2 


1 


1 fii 


n 2 k 2 

u 2 K 


+ n 2 uj c — 


UJ 2 


£02 


$01 ^2 


ulkl 


+ Q.\uj c — 


UJ 2 



(37) 



Finally, to complete the set of pump quantities, we find 
that the velocity in the y-direction is given by 



Ayl _ (Ujk 2 - J^h - (U 2 k 2 - L0 2 )Q 
A x i Sm^oiUjkl + fi 2 ^ c - wg) 



(38) 



with A y 2 given similarly. This is quite different to the 
single ion species case, where A y = 0. 

The behavior of \A XtV \/b is dramatically altered by the 
choice of branch from the dispersion relation; see Fig. 21 
where for brevity we plot A x i, y i only. We see that as 
fco — > 00, the ^-component approaches a finite positive 
value while the y-component approaches zero for the 
upper branch of the dispersion relation. On the other 
hand the x-component approaches large values and the 
y-component passes through zero for the lower branch of 
the dispersion relation. 

Note that the two branches stemming from the up- 
per branch of the dispersion relation become singular as 
fco — > 0. This motivates the need for a normalization 
condition relating the pump field amplitudes of the ve- 
locities A x \^2, Ayi^, magnetic field b, and density i?i,2- 
A reasonable way to proceed is to assume that as we vary 
<x>o the energy density Ed in the pump system is fixed. 
The appropriate relation in a fluid description for general 
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0.1 0.2 0.3 

Normalized Pump Wavenumber v A1 k /a i 



FIG. 2: The ratio of normalized a;- velocity amplitude and 
y- velocity amplitude to magnetic field A x \ tV \/bvAi ver- 
sus normalized pump wavenumber VAiko/£li, with 801 = 
1.1, Q 2 /Qi = 0.1, Pi = 0.8, Bi = 1.5 and U2/U1 = 1/2. The 
solid curve is for the upper branch of the dispersion relation, 
the dashed curve is for the lower branch of the dispersion re- 
lation. Note that on the vertical axis we plot the absolute 
value of this ratio, since we are interested here in the ratio of 
amplitudes. Consequently one of the two dashed curves expe- 
riences a discontinuity in its first derivative when it touches 
the horizontal axis. 



magnetoacoustic waves is (e.g. 24] ), 



B 2 

Ed = — 
1 



2^02 (A 



2 

x2 



K2) 



(39) 



where B — B b, and the summation in the final term 
is over each plasma species. In our approximations the 
electric energy term may be neglected compared to the 
magnetic energy. By switching to dimensionless units 
and setting [IqEd/B 2 = 1, then all our variables are 
given uniquely in terms of the basic plasma parameters. 
By using the expressions for Axi,y\ and A x2 , V 2 we may 
solve for b. The result is 



b = 



V2 



1 1 "". (1+F '. ) "". (1 + F '. ) 

1 1 1.2 2 T2 - 1:2 2 T2" 

"o"Al A l ll V A2 J1 2 



(40) 



ft 



ih. 

-i-! 



B1S0 



[using Fi ; 2 from Eq. (|38|l ] . A plot of b versus wavenumber 
is given in Fig.[3]for each branch u>o(fco)- It is important 
to note that b = 0(yAiko/£li) as fco — > for the upper 
branch. It follows that all the velocity amplitudes are 
well behaved near fco ~ 0. 

The remaining pump quantities to solve for are the 
charge ratios S^l. These may be obtained by imposing 
the charge neutrality condition on both the background 
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0.1 0.2 0.3 

Normalized Pump Wavenumber v k^£2 
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FIG. 3: The magnetic field amplitude b versus normalized 
pump wavenumber VAiko/Q,i due to the normalization con- 
dition JDJ, with 5oi = 1.1, fi 2 /fii = 0.1, fa = 0.8, Bi = 1.5 
and U2/U1 = 1/2. The solid curve is for the upper branch 
of the dispersion relation, the dashed curve is for the lower 
branch of the dispersion relation. 



plasma and on the perturbations to the equilibrium. This 
is justified since we are interested only with frequencies 
that are much less than that of the electron plasma fre- 
quency. This means the electrons have sufficient time to 
respond to the perturbations from the equilibrium and 
neutralize the plasma at each point in space. Thus, on 
the perturbed plasma we impose the condition, 



en e + Z\eri\ + Z 2 en 2 = 0, 



(41) 



where we have ignored the effects of dust (or ion) charg- 
ing, that is, Z 2 is constant. It is then straightforward 
to obtain the charge ratios in terms of the densities, and 
subsequently in terms of the velocity amplitudes and b, 
as follows: 



tfi = <5oi ( 1 + s(b -A x i) cos(fc x) cos(woi) 

LU 



k 



4 0) = ^02 ( 1 + s(b - —A x2 ) cos(fc x) cos(w t) ) • (42) 



This implies that the amount of charge that resides on 
each species is temporally coupled to the other. 

Using the derived set of quantities characterizing the 
pump wave we may, in principle, investigate the possibil- 
ity of excited modes propagating at arbitrary angles with 
respect to the z-axis, as in Ref. However, by remov- 
ing the ^-dependence in the pump waves our calculations 
become much more tractable, with the requirement that 
our excited waves are restricted to parallel propagation 
along the z-axis, see Ref. @- 

For an infinitely extended medium, the perpendicular 
wavenumber of the excited waves is then zero. How- 



ever, this would clearly be inconsistent with wavenum- 
ber matching of a propagating pump with finite perpen- 
dicular wavenumber. We therefore envisage this anal- 
ysis to apply to a standing wave (in the transverse di- 
rection) pump in a finite geometry system, such as a 
cylindrical or toroidal laboratory plasma, or an astro- 
physical magnetic flux tube; by seeking localized solu- 
tions for distances transverse to the magnetic field such 
that k$x = 0(e) in a planar plasma, or kgr = O(e) in a 
cylindrical plasma, we can recover the initially proposed 
form for the magnetic field Eq. (TSJ. In other words, we 
are taking x small rather than ko small, to ensure that 
the fields will look approximately uniform in the perpen- 
dicular direction, i.e. that the wavelength is effectively 
very long on the length scale of interest. In the parallel 
direction, however, the pump fields are uniform, we place 
no such restriction on distances and excited wavelengths, 
and we can apply wavenumber matching of the excited 
wave fields in that direction, as discussed in Section VI. 
Then in the region near the z-axis we have 



B (o) 



o (0) 
P1.2 



(«) 



kox^e 
knx^e 



= B (l + sb cos(to>o£))z, 



kox^e 



(43) 
(44) 



An,2(l + £ — A x i i2 cos(Lo t)), (45) 

UJQ 

£01,2(1 + e(b - ^A xh2 ) cosfwot)), (46) 

UJQ 



correct to O(e). Since the magnetic field has now sim- 
ply a periodic variation in time it resembles a canonical 
parametric pump. The problem then resembles a para- 
metric amplifier, such as a harmonic oscillator with a 
time- varying spring constant. For use in the excited wave 
equations we must also compute the derivatives near the 
z-axis, given here: 



flB(°) 



kftXr^S 




dx 



- kox^e 



koxr^e 



- kQX^e 



0, (47) 
eko{A x \ sin(w £), A yl cos(woi), 0), (48) 
ek (A x2 sm(w t), A y2 cos(oj t), 0), (49) 
= 0. (50) 



am 



dx 



Note that the velocity has a finite derivative, even though 
the velocity itself is zero in this regime. 



IV. PERTURBATION METHOD 

Next we wish to test the stability of the self consistent 
linear solution given in Eqs. I|43|l - (|50|l . now regarded as 
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a finite-amplitude pump wave, to the excitation of waves 
propagating along the magnetic field direction. We must 
however be careful since this is only a linearized pump 
solution. To proceed we will adopt the basic methodology 
of Refs. 0,0]' but shall attempt to refine the argument. 

The basic technique is to perturb each quantity, i.e., 
B(°), v(°), p(°\ fl'°), by an arbitrarily small amount as 
measured by an expansion parameter e' . For example, 
take an arbitrary field, say X, with a known linear rep- 
resentation of the pump field: X^ = X + EX. Then 
Refs. 0, proceed by perturbing this quantity in the 
following way: 



X = X^+e'X'. 



(51) 



However, since X^ has neglected terms of 0(e 2 ), the 
term of 0(s') is even more so negligible in this expansion 
(recall that e' is arbitrarily small). The problem lies in 
the fact that we are ultimately interested in the stability 
of the exact pump solution, say Y"(°) = Ao + eY, with 

Y = X + 0{e). (52) 

Thus, we shall perform the following expansion instead: 

X = Y^+e'X'. (53) 

Upon substitution into our set of nonlinear differential 
equations l|ll|) - (|15[) . we obtain a set of equations of the 
following structure (no approximations): 



sA(Y) + e 2 B(Y) 



e'C(X') + e' D(X') 
Ee'E{Y,X'), (54) 



where A, B,C, D, E are differential operators acting on 
their respective arguments which we may treat as func- 
tions (for brevity we have suppressed their dependence 
on fields other than X). Now eA{Y) + e 2 B(Y) = by 
definition of an exact pump solution. Also, the term of 
0(e' 2 ) can certainly be neglected. Next, we Taylor ex- 
pand E(Y,X') and write 



E(Y,X')=E(X,X') + 0(e) 



(55) 



We are then permitted to neglect this term of O(e) when 
inserted back into Eq. I)54[l. This gives us a set of equa- 
tions for the perturbation X' , in terms of the X, to suf- 
ficient accuracy, i.e., 



C(X') =eE{X,X') 



(56) 



Note that X has served two purposes. First, it gives us 
an approximate description of the corresponding exact 
solution = Xq + eY, and second, it allows us to 

investigate its stability without ever having to find Y . 



V. EQUATIONS OF MOTION AND NATURAL 
MODES 



Our task then is to ascertain the primed variables 
which are the excited fields. As stated previously, we 



are interested only in plane waves that travel parallel to 
the z-axis. In this case, any longitudinal components will 
decouple completely and merely describe a linear acous- 
tic wave, so we may set the z-componcnt of the primed 
velocity to zero. Hence, both primed velocities have the 
following form: 



'1,2 



Kl,20M)> u yl,2fo*))O)- 



(57) 



Thus, V • = V • v' 2 = 0, and from each of the two 
continuity equations we have 



dp'1,2 
dt 



-ep' 12 [V • vi, 2 ] 



(58) 



where we have imposed the kox = O(e) [or k$r = 0(e)] 
condition. Hence we have p[ = p' 2 = 0, i.e., there are no 
perturbations in the densities of the excited fields. From 
this it also follows that 6[ = S' 2 = 0. 

By taking the derived expressions for B, v, p and 8 and 
substituting into Eqs. (|13fl & i|14|) . and using the stability 
analysis procedure conveyed in Eq. I|56|l . we obtain 



dt VxKxBoH^Vx ^ 



(59) 

kox^e 



<9v' 1 

Poi $oi-^r ~ P02^2(vi - v 2 ) x z (V x B') x B 

at po 



-5 m p m v' xl — + n 2 p 2 (vi - v 2 ) x z 
ox 



(60) 



knx^e 



Here the terms of O(e), which occur on the right hand 
side, should be thought of as driving terms from the 
pump wave. Note that there exist two further equations 
of motion under the index interchange 1^2. Note also 
that there are no acoustic terms on the left hand sides of 
these equations because of the decoupling of the longitu- 
dinal motions: however there are still acoustic influences 
in the pump fields on the right hand sides of the equa- 
tions. 

Now in order to treat the x and y components of the 
B and v vectors on an equal footing we form a complex 
vector out of each component of these partial differential 
equations, utilizing the variables v± = v x ±iv y . The spa- 
tial variation is assumed to be periodic in the z-direction 
with wavenumber fc, as follows: exp(ifcz). Given this, the 
following pair of linear ordinary differential equations are 
obtained: 

dB± kB dv±i 

- ikB v ±1 T 



dt 



xe 



Bobkv±i + k 



fii dt 

BokoAyi , , , 

v xl ± k A yl B x ) cos(woi) 



T A xl (kaB' y + ik ~^-^j sin(o> t) 

OoiPoi^r— ± «P02"2(w±i -v±2)- i B± 

dt p 



(61) 
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5oiPoikov' xl (A xl sin(w t) ± iA yl cos(u t)) 



T i^2P02 — A x2 cos(uj t)(v±i - v±s) 



(62) 



Let us define the sense of polarization in reference to the 
screw sense of the fields in the direction of propagation in 
the z-direction. Then t>+ = v' x + iv' y corresponds to a left 
hand circularly polarized wave for positive frequencies 
and a right hand circularly polarized wave for negative 
frequencies, while V- — v' x — iv' corresponds to a right 
hand circularly polarized wave for positive frequencies 
and a left hand circularly polarized wave for negative 
frequencies. Also, B± is defined similarly. With these 
definitions v x ,v', B' x , B' y may be eliminated entirely from 
these expressions. Note that the presence of a second 
moving ion species introduces coupling through terms of 
the form: v±\— v± 2 , and an additional term on the right 
hand side in Eq. (|62J) due to density variations. 



In an analogous way to the classic driven pendulum 
problem [25| . we have found that the velocities of the 
two species satisfy a pair of generalized Mathieu equa- 
tions. Furthermore, the two species and magnetic field 
are strongly coupled. To proceed, we move from the time 
domain to the frequency domain, under the Fourier trans- 
form: 



v ±1 {t)e luJt dt, 



i r°° 

V ±2 (co) = -= / v ±2 (t)e^dt, 

V 2-7T J-aa 



(63) 



where u> is the frequency of the excited waves. Using the 
linearity of our pair of differential equations, we compute 
the Fourier transform of both. Upon eliminating B± and 
V± 2 (uj) in favor of V±i(lo), we obtain: 



F±(lo)V ±1 {lo) = e 



where, 

F±(to) 

M ± (A 1( 

N^A 2 ,D u e) 



+ M ± {-A x1 ,+A x1 ,uj + ,-1)V +1 (uj + 
+ M±(+A x1 ,-A x1 ,uj-,-1)V + i(uJ- 
+ M±(-A x1 ,-A x1 ,uj + ,+1)V- 1 (uj + 

+ M±(+A r i,+A E i,w_,+l)U_i(w_ 



+ V±i(w_)) + \Pd!») (Vufa+) + V±&-)) 
+ N±{~A x2 , +A xl ,-l)V +2 {uo+) 
+ N ± {+A x2l -A xl ,-l)V +2 (tJ-) 
+ N ± (-A x2 , -A xl , +1)U_ 2 (cj + ) 

+ N±(+A x2 ,+A xl , +l)U_ 2 (w_) 



Soilo 2 +v 2 A1 k 2 (It — 
koV x2 



fliSoi 

± — : U) 



Q 1 
b 



± — : U> 



1 - 



it 



'112 



1 =P — 

1 t no 



'02 



CJ 



it 



fco 

4 



(5 0lW ± ^-)(A, ± A yX ) ± (-5 01 £ + e(l - <S„i)fii) ( 1 - 



1 



fc(lT^) 



±-e(l - 5 i)fii(-Di T A yl ) - — ^ — , , 



(64) 

(65) 

(66) 

(67) 
(68) 



and w± := cj ± Wo- Eq. (|64|l is obviously not an algebraic 
expression from which we can uniquely obtain the dis- 
persion equation between w and k. Instead it provides us 
with a functional relationship between the Fourier trans- 
forms with arguments to, lu ~ luq and to + luq. 

Nevertheless, we may proceed by noting that the right 
hand is 0(e), which is small. For the left hand side to be 
small we must ensure that the frequency is near a root 
of the polynomial on this side with a correction of 0(e), 
i.e., 

to = X + 0(e), (69) 



where \ is defined to satisfy, 

Fdx) = 0- (70) 

These define the naturally occurring modes in the ab- 
sence of pumping, (i.e., if e = then to = x)- The solu- 
tions of these algebraic expressions give the dispersion re- 
lations characterizing plane transverse waves propagating 
parallel to the background magnetic field in a two species 
plasma. The last term of F± is only present when a sec- 
ondary moving ion-like species resides in the plasma. It 
is an additional Hall-type term, associated with the rela- 
tive motion of the two charged species. In the one species 
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limit this relation simplifies to 



X 2 + v 2 A1 k 2 



0. 



(71) 



and we are left with the familiar fast and slow (ion- 
cyclotron) transverse waves 

The important feature of the two-species result in 
Eq. H7Ufl is that it is a cubic in \ an d hence another 
mode of excitation has been added. Moreover, since we 
may choose F + or F_ it follows that there are a total of 
six solutions; three are left hand circularly polarized and 
three are right hand circularly polarized. However, only 
three of these are physically different. This is because 
if x satisfies F + (x) = 0, then — x satisfies F-(— x) = 0, 
and vice versa. We may then concentrate on the positive 
frequency solutions. 

If the secondary species is positively charged then there 
are two left hand modes of excitation 01 ■ These we de- 
note by u>L a and wl/, where the "s" denotes slow, and the 
"/' denotes fast. There is also a single right hand mode, 
denoted by lor. The corresponding dispersion relations 
are plotted in Fig.^Ja). Note that the right hand mode 
intersects with the fast left hand mode. Also, the two 
left hand modes both experience a resonance as k — > oo. 
For LULf it is given by f2i, and for ujl s it is given by O2 
(where we have assumed Ox > ^2). 

For a negatively charged secondary species there are 
also three modes of excitation, however the combination 
of polarizations has changed. In this case there are two 
right hand modes u>Rf and ujr s , and a single left hand 
mode u>l', see Fig.0Jb). Here ojr s has a resonance at 
— 0,2, and wl has a resonance at Q,\. In the dusty plasma 
case, further properties of these modes were investigated 
inRefs. HISS- 

Note that the upper curves (fast modes), regardless of 
the sign of the species, experience a nonzero cutoff fre- 
quency which coincides with that of the pump frequency, 
see Eq. (jUJ. 



VI. GROWTH RATES OF PARAMETRIC 
INTERACTION 

From the above analysis we see that x may be any 
one oi ujls, u>Lf, for a positively charged secondary 
species, or ujr Si ^>rj, for a negatively charged sec- 
ondary species. Since the frequency describing the ex- 
cited wave is a perturbation of a natural mode frequency 
we denote the change by <j) which is allowed to be com- 
plex, where 



CO 



(72) 



Now by returning to the Fourier transform relationship 
in Eq. I|64ll we can see what effect <p (i.e., pumping) will 
have on V±. Without any pumping (i.e., s — 0) the spa- 
tial solutions would be a linear superposition of three 
pure monochromatic exponentials. Hence the Fourier 
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FIG. 4: The normalized excited wave frequencies x/^i ver ~ 
sus normalized wavenumber VAik/Sli of the three natural 
modes, (a) The second species is positively charged, with 
<5oi = 1.1 and = 0.1. (b) The second species is nega- 

tively charged, with 801 = 0.9 and fb/fii = —0.1. 



transform would be a sum of three Dirac-(5 functions; see 
Fig-EIa) for a representation. The effect of pumping is to 
modify this and provide frequency shifts. The modified 
solution will have some broadened frequency spectrum, 
since it is the case that Dirac-(5 functions do not solve 
the Fourier transformed equations in the presence of a 
pump. These types of modifications will occur near each 
natural frequency; see Fig. 5(b) for a simplified represen- 
tation of this effect. In fact an exact solution would be 
more complicated with the identification of harmonics 
etc, but this figure should illuminate the key feature of 
the interaction. 

This is described as a parametric interaction between 
the excited fields due to the pump fields. In order to 
proceed we consider the case in which the interaction is 



10 



> 




s 
> 
E 
o 

M— 

w 
as 



-to L 











Excited Frequency to 




(b) 








) 






1 



(o 

u Rs 

Excited Frequency m 



FIG. 5: A representation of the Fourier transform V(ui) versus 
frequency lo. (a) No pump, giving three Dirac-delta functions, 
(b) The effect of a pump, broadening the spectrum. Note 
that we have chosen to illustrate the idea using only the right 
hand modes (i.e., V-), for a plasma with a negatively charged 
secondary species. 



greatest, as in Refs. 0,0- It is clear that this effect will 
be greatest when the right hand side of Eq. l|fjl|) is large, 
which is the resonance condition. This occurs when one 
of the Fourier transforms on the right hand side is reso- 
nant, in other words when their arguments are near one 
of the x' s j as conveyed in Fig.0 Since the right hand 
side of Eq. I|t)4[) involves terms of the form: V(lo+) and 
V(lu-), it follows that there will be a large parametric 
interaction when lo + loq or lo — loq satisfies this. Now, as 
mentioned earlier we are restricting our attention to pos- 
itive lo (without loss of generality). Under this condition, 
it is found that only lo — loq can satisfy this. Moreover, 
there are three choices for this interaction: a left hand 
mode interacting with a right hand mode, a left hand 
mode interacting with a left hand mode, or a right hand 



mode interacting with a right hand mode. 

To illustrate, suppose the secondary species is nega- 
tively charged and ui « ljl (i.e., \ = lol)- Then we 
can have lo — loq w — lOr s or lo — luq w —lorj (left hand 
mode interacting with right hand mode) or lo — loq ps — lot, 
(left hand mode interacting with left hand mode). In 
other words, the resonance condition is that the pump 
frequency satisfies 

LOq = LO L + LO Rs , LO = LO L + LO R f, or LO = 2lo Li (73) 

respectively, which is a statement of conservation of en- 
ergy. That is, the resonance condition describes the cou- 
pling of a pump wave with two daughter waves. 

Moreover, if lo w lor s then we can have lo — loq ps — lol, 
lo—loq w —<-^Rf or lo— loq w — lor s and the same conditions 
for resonance apply. Similar rules apply for lo w lor f and 
for a plasma where the secondary species is positive. 

Let us now address the issue of the treatment of the 
wavenumber k in our pair of interacting waves. First, 
note that in the z-direction the pump wavenumber fco is 
zero. For the natural modes we have that a wave given 
by (k, — x) is physically equivalent to a wave given by 
(— fc, x)- Thus our two interacting waves will be given by 
(k,xi) and (— k,X2), where Xi,X2 are any combination 
of LOL,LOL S ,LOLf ,LOR,LOR S ,LOjif . Hence the wavenumbers 
associated with each interacting wave are equal and op- 
posite, i.e., 



fci + k 2 = 0, 



(74) 



which is a statement of conservation of momentum in the 
z-direction. The fact that the right hand side is zero re- 
flects the spatial uniformity of the pump wave. This tells 
us that if the approximately spatially uniform standing 
pump wave decays, then it does so into two daughter 
waves of equal wavelength travelling in opposite direc- 
tions. Since the wavenumber magnitudes are equal we 
may just consider the wavenumber k without referring 
to the sense of polarization. This is a direct generaliza- 
tion of the previous investigations in 0, H, El • 

With the above framework we now proceed to solve for 
<f> from the Fourier transform relationship. We discuss 
the method involved in obtaining (j) in the case where 
a left handed wave interacts with a right handed wave, 
denoted (f>LFt- Here we may assume, without loss of gen- 
erality, that x = Xl in Eq. 172fl . where xl is any one of 
LOL,LOL S ,LOLf . This will interact with — xb. (i- c -, X ~~ 
may be any one of — lor Si —Lo R f, —lor). 

Following the methodology of Ref. [28| we form another 
equation corresponding to (|64fl under the transformation 
lo — > lo — loq. In the resulting two equations, we neglect 
the obviously nonresonant terms V(lo + ) and V{lo — 2oj ). 
Additionally, due to our particular choice of interaction, 
i.e., left-right, we neglect V+{loJ) and V-{w), and retain 
only V+(u>) and V+(w_). We find the following: 



F+(xl + £<I>lr)F-(xr - e(j) LR ) 
= e 2 [Y+{-xr) + J+(+A xl ,+A x2 , 
x [Y+{+ X l) + J+(-A xl ,-A x2 , 



-A xU 
-A xl , 



-Xn, 

'XL, 



-1)] 

•1)] ,(75) 
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where 

J ± (A 1 ,A 2 ,D 1 ,C,e) 




,(76) 



IT-* 



^(^D^e) 



+M ± (A 1 ,£» 1 ,e,e). (77) 



In these expressions we implicitly have that £ and w, are 
related by: £ — > xl => w — > ~xr and £ — > X-R => w — > 

-XL- 
Ill the one ion species analysis, it is at this point that 
the conservation of energy and momentum rules are used 
to explicitly obtain xl, XR, k. The procedure is to use F± 
to solve for Xl,r m terms of k. These solutions are then 
added together and equated to u)q. The resultant expres- 
sion is then solved for k and subsequently substituted into 
the expressions for Xl,r- This process explicitly shows 
that the resonance condition, for a given value of pump 
frequency wq, uniquely determines the frequencies and 
wavenumbers of the excited waves. Now, this is all pos- 
sible since the solution of F± in the one species limit is 
merely a quadratic, see Eq. I|71[) . However, in our case 
we must solve a cubic. As such, the resulting expressions 
for xl, X-R, k are too complicated to be reproduced here. 
Let us just note that it is still true that loq determines 
them uniquely. 

In order to solve for 4>lr in Eq. I|75|) we Taylor expand 
F + (xl + £4>lr) around xl to order e, then use the fact 
that F + (xl) = 0, to obtain 



F+(XL 



£<Plr 



^lrH( X l) 



(78) 



where 



0,1 



M«2-x) 2 



(79) 

while F_(xr~ £ 4>lr) can be treated similarly. Hence 4>lr 
is obtained from: 

4>Ir = [Y+(-Xr) + J+(+Axi, +A x2 ,-A x1 ,- X r, -1)] 
x [Y+(+ X l) + J+(-A xl ,-A x2 , +A xl , +xl, -1)] 
x [H(xl) HI-xrT 1 ■ (80) 

Although this was derived assuming we were near a 
natural frequency of a left hand mode, we obtain the 
same expression for 4> 2 for the right hand case (viz. <j> RIj ) 
with the modification that 4>rl is of the opposite sign. 

This technique may also be applied to both the left- 
left interactions and the right-right interactions. This 
includes the possibility of a slow wave interacting with 
a fast wave (e.g., u>l s and Wl/), and waves of equal fre- 
quency interacting. These sorts of interactions occur only 
because of the velocity derivative terms in (|48|l & l|49|) . In 
obtaining cf> we must retain and neglect the appropriate 
selection of Fourier transforms according to the choice 



made. Here we introduce xli and XL2 as two (possible 
equal) left hand modes, and x_ri and xri as two (possible 
equal) right hand modes. The corresponding frequency 
changes are obtained from 



f L L = [J+{+Axl,+A x 2,+A xU -XL2,-l)} 

x [J-(-A x i,-A x2 ,+A x i,+XLi,-l)} 

X [H( X Ll) H(- X L2)V\ 
f>RR = [J-{+A xl ,+A x2 ,-A xl ,-XR2,-l)] 

x [J+(-A x1 ,-A x2 ,-A x i,+xri,+1)} 



(81) 



[H(xri) H(- X r 2 )Y 



(82) 

for the left-left and right-right interactions, respectively. 

It is worth noting that in other types of parametric 
problems, such as modulational and beam instabilities, 
a graphical approach is useful in classifying the instabil- 
ities; they occur where two normal mode lines cross and 
reconnect [2t| H(j • In those cases the number of interact- 
ing modes is finite, and the nonlinear dispersion relation 
can be obtained in closed form. In our case, however, we 
cannot obtain the nonlinear dispersion relation in closed 
form, due to the infinite number of interacting modes; 
the decay interaction of the two excited modes with the 
pump is postulated in the first approximation, and the 
classification of the instability is relatively straightfor- 
ward. 



VII. NUMERICAL ANALYSIS 

From the expressions for <f> LR , </>ll, ^rr it is easy to 
see that these quantities are real valued when we are in 
a frequency-wavenumber regime such that a pump wave 
can propagate. Hence is either purely real or purely 
imaginary. If <f> is real then what we have found is a 
correction to the natural mode frequency which adds an 
extra Fourier component to the field expansions. If <p is 
imaginary then the daughter waves will either undergo 
exponential growth or decay. Moreover, since there will 
exist both a positive and a negative imaginary frequency 
solution and since the frequency is conserved, this implies 
that one of the daughter waves will grow and the other 
will decay. In such cases the presence of one exponentially 
growing daughter wave implies that anywhere in which 
4> 2 < is a region of instability. In this case 4> is the 
growth rate. 

In Fig. we plot the full complement of (normalized) 
squared frequency changes. That is, we plot (4>lr/Qi) 2 
versus pump frequency wo/^i for the upper "{/" and 
lower "L" branches of the pump dispersion relation for 
each mode of interaction. This figure is for a warm 
plasma in which the second heavy species is positively 
charged. The same basic features are present in the neg- 
atively charged case. All curves labelled "s", i.e., "Us" 
and "Ls" , correspond to the combination between ujr 
and ujls, that is the slow interaction. The curves labelled 
"/", i.e., "[//" and "£/", are for the fast interaction be- 
tween lor and u>l / ■ 
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FIG. 6: The squared normalized frequency change (4>lr/&i) 2 
for the left-right interaction plotted against pump frequency 
too/fti- In both plots the second species is positive, with 
fa = 0.8, B x = 1.5, U 2 /U 1 = 1/2 and Q 2 /Q 1 = 0.1. The 
number densities vary in (a)-(b), with 601 = 1.02 in (a), and 
Sqi = 1.002 in (b). Also, "Us, Uf, Ls, Lf" denotes the upper- 
slow, lower-slow, upper-fast and lower-fast combinations, re- 
spectively. Here "N" denotes the single species result. 



In Fig.|Sfa), where <5rn = 1-02, the full spectrum of per- 
turbations to the natural modes is seen. The U f interac- 
tion is unstable and monotonically decreases, approach- 
ing the single species result N for large ujq. Hence this is 
indeed consistent with the single ion analysis, given by 



,A 1,4 



1. 



(83) 



The Lf interaction is also unstable, however it experi- 
ences a minimum (corresponding to a maximum growth 
rate), about which it turns over and approaches zero. For 
our choice of parameters, in this figure, this occurs near 
uiq ~ 0.2 ili. Moreover the intersection with N occurs at 



FIG. 7: The squared normalized frequency change (^ll/Qi) 
for the left-left interaction plotted against pump frequency 
loo/SIi. Here the second species is positive, with 0i = 
0.8, Bi = 1.5, U2/U1 = 1/2, fia/fti = 0.1 and 5 i = 1.002. 
Also, Uss, Lss, Uff, L f f, Usf, Lsf denotes the upper-slow- 
slow, lower-slow-slow, upper-fast-fast, lower-fast-fast, upper- 
slow-fast and lower-slow-fast combinations, respectively. 



u>o ~ 0.45 f2i. Next, we see that the Us interaction, is 
small and positive, and hence it is weakly stable. Note, 
however, that if higher order terms are included in the 
field expansions then this may in fact prove to be weakly 
unstable. Finally, we see that the Ls interaction, which 
is only present in a warm plasma, is strongly unstable, 
with a growth rate that approaches infinity at a rate con- 
siderably faster than N or Uf. 

As we decrease the number density of the second 
species in Fig.l^b) to <5ni = 1.002 some interesting fea- 
tures appear. The Uf curve shows a sharply varying 
growth rate at just above the cutoff frequency. The Lf 
minimum shifts to a lower value (lying beyond the range 
of the plot). The Us interaction has actually become 
unstable in the frequency region just above the cutoff 
frequency. Also, the manner in which our bi-ion result 
approaches the single species result as Sqi — > 1 is inter- 
esting: the envelope of the three modes Ls, Us, and Uf 



< 



< 



LJ 



5- 2tu r and 



form in the domains ujq 
loo 2 lo c , respectively. 

In Fig.[7| we plot the range of left-left interactions for 
the positively charged secondary species case. As can be 
seen in Fig.^a) this case has two left handed branches 
(slow and fast). Hence we may form the slow-slow, fast- 
fast and slow-fast combinations. This combined with the 
lower and upper branches of the pump dispersion relation 
gives a total of six interactions. Of these six we find that 
all are unstable, except the "is/" combination which 
experiences a point of maximum stability at just below 
2ut c . Note that the "Lss" interaction decreases without 
bound as uiq — > 2f2i This is acceptable given that for the 
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VIII. CONCLUSIONS 
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FIG. 8: The squared normalized frequency change (4>rr/Qi) 2 
for the right-right interaction plotted against pump fre- 
quency ujq/Q,\. Here the second species is negative, with 
Pi = 0.8, Bi = 1.5, U2/U1 = 1/2, fia/fti = -0.1 and 
5oi = 0.998. Also, u Uss, Lss, Uff, Lff, Usf, Lsf" de- 
notes the upper-slow-slow, lower-slow-slow, upper-fast-fast, 
lower-fast-fast, upper-slow-fast and lower-slow-fast combina- 
tions, respectively. 



slow-slow case the excited wavenumber k experiences a 
resonance as u> — > fi\ (with pump wavenumber small). 
Although it lies below the range of the plot the "Lff" 
interaction starts at a finite value, when ujq = 2uj c . Also 
note the curious feature wherein the U U sf" interaction 
has a corresponding growth rate which is large in the 
range to c < u>o ^ 2 lo c only. 

Finally, in Fig.|H| we plot the range of right-right in- 
teractions for the negatively charged secondary species 
case. Again by inspection of Fig.^b) it makes sense 
that there are six values of 4> 2 RR - This plot shares many 
features in common with the previous plot. However, an 
important change is that the Lss interactions has a min- 
imum and approaches zero as k — + 00. Also, the Usf 
maximal growth rate has shifted to much higher values 
of pump frequency. The Lsf interaction has stability as 
luq increases. 

Hence, we have seen that the introduction of a second 
species adds a range of extra behaviors, even with a small 
number density. Note that in the special case where the 
frequencies of different modes coincide, the selection of 
resonant Fourier transforms to generate these plots is ac- 
tually invalid. In particular, this occurs at the point seen 
in Fig.QJb) where ui R = uJzf- At that point there are 
actually more terms resonant than have been accounted 
for. However, since this occurs on a set of measure zero, 
it has been ignored. 



This paper has outlined an investigation into a cer- 
tain class of instabilities in a bi-ion or dusty plasma, al- 
lowing full mobility of all species and the inclusion of 
pressure. The results should be applicable to a wide 
range of laboratory and space plasmas where secondary 
ion species or charged dust grains occur. The immediate 
consequence of the presence of an extra heavy species is 
its increase or reduction of the number of free electrons 
in the plasma through the charge neutrality condition. 
We obtained a pump wave with spatial variation trans- 
verse to the background magnetic field, which contained 
an upper and lower branch, and made the approxima- 
tion of spatial uniformity of this pump. The stability of 
this structure was investigated by perturbing the pump 
wave, as we looked for parametric instabilities wherein 
the pump wave generates excited waves. Three natural 
transverse modes were found to be excited propagating 
parallel to the background magnetic field; two left handed 
when the second species is positive and two right handed 
when the second species is negative. Two out of each of 
these three are the modified fast and slow (ion-cyclotron) 
modes. There were three basic types of interactions al- 
lowed; left-right, left-left and right-right, each of which 
has several combinations dependent on both the choice 
of branch from the pump wave (upper and lower) and 
the choice of natural mode (slow and fast), giving rise to 
twelve interactions in total for any particular plasma. 

We derived the frequency changes to the natural 
modes, corresponding to growth rates when the perturba- 
tion to the frequency was imaginary. We can summarize 
the main features of the instabilities for the left -right 
interactions (opposite sense circular polarization) as fol- 
lows: 

(i) Of the four combinations allowed, three are unstable 
when the number density of the second massive species 
is substantial, and all four are unstable when the number 
density of the second massive species is low. 

(ii) As the single species case is approached (Si — ► 1), 
the envelope of three different growth rates approaches it 
in the domains uiq < u c , uj c < cl>o < 2 uj c and loq > 2u> c . 

(iii) The two interactions corresponding to the lower 
branch of the pump dispersion relation are strongly un- 
stable, with the Lf experiencing a local maximum and 
Ls unbounded. 

Next we investigated the left-left and right-right inter- 
actions (same sense circular polarization): 

(i) These interactions involve six and two combinations 
respectively when the second species is positive, and the 
reverse when the second species is negative. This fol- 
lows simply from the assortment of handedness of natural 
modes. 

(ii) Of these only the Lsf interaction is stable, while 
the others display interesting behaviors. In particular, 
in the positive case, the Lss curve approaches — - 00 as 
luq — > 2lu c where the excited wavenumber k experiences 
a resonance, while its growth rate is maximal at ujq just 
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below 2\u> c \ in the negative case. 

(iii) In contrast to the left-right interactions, the upper 
pump mode can give the strongest instabilities. The Usf 
interaction has a large growth rate between oj c and 2oj c 
in the positive second species case, and at higher values 
of pump frequency in the negative second species case. 

In several of the above interactions the growth rate 
starts at a nonzero value. This occurs in the cases where 
the starting value of ojq have a corresponding nonzero 
value for k , so the velocity derivative terms are nonzero. 



A possibility for further work would be to allow the ex- 
cited waves to be fully oblique, rather than simply par- 
allel to the background magnetic field. In doing so we 
would be able to remove the long pump wavelength as- 
sumption, which may lead to interestin g fe atures. In the 
single ion case this was treated in Ref. [TlJ , albeit with- 
out pressure. It is anticipated that implementation of the 
accompanying additional wavenumber matching rules, in 
both the z and x directions, would be algebraically com- 
plicated. 
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